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Abstract 

A method for determining the dependent K*^ spin amplitudes of K*^pi~ 
decays through a maximum likelihood fit to data is presented. While current 
experimental techniques extract a limited set of observables in bins of cp', our 
approach allows for the determination of all observable quantities as continuous 
distributions in <p. By doing this, the method eliminates the need to correct theory 
predictions of these observables for averaging effects, thus increasing the sensitivity 
to the effects of physics beyond the Standard Model. Accounting for the symmetries 
of the angular distribution and using a three parameter ansatz for the (p dependence 
of the amplitudes, the precision of the angular observables and the sensitivity to new 
physics is estimated using simulated events. These studies are based on the sample 
sizes collected by the LHCb experiment during Run-I and expected for Run-II. 
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1 Introduction 


Rare b —)■ processes are suppressed in the Standard Model (SM) as they can proceed 

only via electroweak penguin or box type diagrams. As-yet undiscovered particles could 
give additional contributions with comparable amplitudes to those of the SM processes, 
and such decays are therefore sensitive probes of new phenomena. 

The angular distribution of the system in —)■ decays is of 

particular interest, as it can be described by a number of measurable quantities which 
are sensitive to new physics and can be precisely predicted in a given physics model. 
Theoretical predictions for such observables are particularly precise in the range of dimuon 


invariant mass squared, 1 < < 6 GeVyc"^ The potential of the )■ 

decay as a probe of New Physics (NP) has resulted in numerous proposals for observables 

with varying levels of theoretical precision, as described in Refs. 

The dominant uncertainty in the predictions of —)■ observables, is at¬ 

tributed to the calculation of the B^ —)■ K*^ hadronic form-factors. The LHCb collaboration 
have determined a number of observables which are designed to have a reduced depen¬ 
dence on these form-factors 15,^. The measurement of the observable 13 , which is 
designed to ensure the cancellation of the hadronic B^ —)■ K*^ soft form-factors at leading 
order, exhibits a local tension at the level of 3.7 a with respect to the SM prediction of 
Ref 13 . This measurement has been interpreted as an indication of a new, heavy, vector 
particle [6,17-29 or as a consequence of previously unaccounted for QCD effects 30 -32 . 
In addition to potentially unexpected QCD effects, there are several other factors that 
limit the sensitivity of the data to new physics effects. As detailed below, these include the 
omission of certain symmetry relations between the observables in experimental analyses, 
the interference between P- and S-waves of the K~7t^ system, and the binning in g^. 

The LHCb collaboration have published a number of separate analyses of B^ — 


decays to determine different observables 15,33 . Recently, the LHCb collabo¬ 
ration also presented a measurement of all observables related to the K~tt^ system in a 
P-wave state [T^. In all these measurements, the simple relations between the P-wave 


observables, arising due to the symmetries of the angular distribution, are implemented. 
However, the remaining complex relation between P-wave observables 


11,34 


IS not 

exploited. This results in redundant parameters being determined, reducing the overall 
experimental precision of the observables. 

One of the dominant systematic uncertainties for the experimental determination of 
B^ —)■ observables is the lack of knowledge of the S-wave contribution to the 

predominately P-wave system. An S-wave contribution can induce a bias on the 

(P-wave) observables and dilute the experimental sensitivity to such observables 35 . Until 


recently, the experimental results included a systematic uncertainty to cover this. In 
the latest LHCb analysis 16 , the S-wave components are explicitly included. However, 


further symmetry relations between the P- and S-wave components are not accounted for, 
resulting in two redundant S-wave parameters 36 . As in the case of the omitted P-wave 


symmetry relations, this redundancy dilutes the experimental precision of the observables 
of the angular distribution. 
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The precision of global fits to existing measnrements of —)■ transitions 

are also limited by the fact that the measnrements are performed in bins in Com¬ 
parison between theoretical predictions and experimental measnrements then reqnires 
the integration of theoretically predicted observables over experimental bins. This 
integration has the effect of dilnting variations and introdncing dependence on other, 
potentially more poorly predicted, observables. For example, the angular terms in the dif¬ 
ferential decay rate that involve the theoretically clean observable have a 
prefactor, where Fl is the longitudinal polarisation fraction of the K*^ and is a poorly 
predicted observable. This results in the experimental measurements being sensitive to, 
J ^FL(g^)(l — Fi,{q^))P^{q^)dq^ . Although the angular distribution is also sensitive to 
/ this does not enable J P^{q^)dq^ to be computed and hence does not allow 

the full exploitation of the cancellation of the form factors at leading order for which P^ 
was designec0 

In this paper, we propose a method of analysing 5° —)■ decays that allows 

the determination of all of the K*^ amplitudes as a parametric function of This 
method allows the formation of any observable from a single fit to the data from a given 
experiment, including the full experimental correlations. The method also allows the 
S-wave related amplitudes to be determined, removing the need for any experimental 
systematic uncertainty from S-wave contamination. The g^ shape information and the 
application of all the symmetry relations of the angular distribution, result in a substantial 
gain in sensitivity to new physics effects and completely remove the g^ averaging problem 
mentioned above. The method can be used with the —)■ sample that is already 

available at the LHCb experiment. 

The paper is organised as follows: Sec. |2]describes the decay rate and angular observables 
of 5°—)- K~Ti^jjpiji~ transitions; Sec. [^describes the method for extracting the K*^ helicity 
amplitudes from a single fit to data; Sec. [^presents the results of the method using simulated 
B^ —)■ decays with sample sizes equivalent to those obtained by LHCb during 

Run-I and projected for Run-II of the LHC; finally Sec. [^compares the sensitivity of the 
amplitude fit to other methods for extracting the K*^ angular observables. 


2 The differential decay rate 


The differential decay rate of the B^ meson, for a Kn system in a P-wave configuration 
and ignoring scalar contributions to the dimuon system, is given by 


^Experimental measurements 15 16 33 also enable the g^-averaged to be computed using the ratio 
{S^)/ yj (El)(1 — (El)), where {S^) and (Fl) denote ( 7 ^-averaged quantities. This ratio is not the same as 
the optimal observable — Fp)). 
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[ Jis sin^ 9k + Jic cos^ 9k + J2s sin^ 9k cos 29^ + 

OZTT 

J 2 c cos^ 9k cos 29i + J 3 sin^ 9k sin^ 9i cos 2<j) + 

J 4 sin 29k sin 29^ cos 0 + J 5 sin 29k sin 9^ cos 0 + ( 1 ) 

Jqs sin^ 9k cos 9i + J 7 sin 29k sin 9£ sin 0 + 

Jg sin 29k sin 29^ sin 0 + Jg sin^ 9k sin^ 9^ sin 20 ]. 

The angular observables, Ji{q^)-, depend on six dependent complex amplitudes, A^’^, 
representing the three polarisation states of the K*^. The conhguration of 
a longitudinally polarised K*^ and time-like polarised dimuon system is suppressed and 
therefore safely neglected. The labels L and R refer to the chirality of the dimuon system. 
The various Ji{q^) observables are given by 


d cos 9^ d cos 9k d 0 dg^ 


Jis = [\Ai\^ + |klf|2 + (L ^ /?)] + ^Re{AiAf* + A^Af*) 

4-777 ^ 

Jie = I4r + Idlo^'r + ^[2Re(4Tf)] 

J2s = + l^f P + (L -)■ R) 

J 2 C = —l3l \Aq\‘^ + (L —)■ i?) 

^3 = f ikiir-ikifr+(L^R) 

r 1 

J 4 = ^ Re(kl^Tf) + {L^R) 

J 5 = V2/3^ Re{A^Ai*) - {L ^ R) 

Je, = 2/3^ Re(RfTi*) - {L ^ R) 

R = V2R lm{A^A^*) -{L^R) 

Jg = ^ + {L^R) 

Jg = 0 ; Im(Rf Ti) + iL^R) 


( 2 ) 
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with = (1 — 4m‘j^/q‘^) and (L — R) denotes the same term bnt with flipped chirality. 
In the limit of S> 4m^ the various Jj coefficients can be trivially related by J 2 c = —Jic 
and J 2 s = </is/3. An additional, more complicated relation exists between J 2 c and the rest 
of the angular observables as noted in Ref. 


11 and explicitly given in Ref. [34| . 

While the differential decay rate in Eq. [I is defined for the decay of the meson, the 
decay of the can be given in complete analogy, by starting from Eq. and performing 
the substitution Jj —)■ Jj, following the angular convention described in Refs. 


15 33 This 


convention allows the sum of the decay rates of B^ and mesons to be written in terms of 
sums of Ji and R angular observables by simply performing the substitutions R —)■ (R + R) 
into Eq. (^. It is therefore convenient to define CP-averaged observables Si, as discussed 
in Ref. (T , 


S^ = 


Ji ~\~ Ji 


(dr + dr) /dg2 ■ 

In a similar way, the CP-asymmetric observables, Aj, are defined as 

Ji Ji 


( 3 ) 


A. = 


(dr + dr) /dg2' 


( 4 ) 


2.1 S-wave interference 


The expression in Eq. Q assumes that the system is in a P-wave configuration, as 

is the case for the A'*°(892) vector meson. To account for a K~tt^ system in an S-wave 
configuration, the spin-amplitudes need to be modihed to account for the presence of the 

L R 

S-wave amplitudes Aq ’ 


^00 


In previous experimental analyses 15,33 the presence of an S-wave contribution was 


accounted for by assigning a systematic uncertainty. The method presented here enables 
the determination of the S-wave amplitudes using a modified version of Eq. 


d^T d^T 

d cos 6^ d cos 9k d(/) dg^ d cos d cos 6k d0 dg^ 

- cos29i) +J”icCos9k{ 1 - cos29i) + 

32vr ( 5 ) 

J '4 sin 29i sin 9k cos (j) -f J '5 sin sin 9k cos 0 -f 
J'Y sin 9£ sin 6 k sin (p + J'^ sin 29£ sin 6 k sin (p ], 

with 
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J'7 = 2\ ^ 
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ReiA^,Ad) + iL^R) 
Re(/l^o^f) + {L^R) 
Re{A^,Ad) -(L^R) 

Im(R^o4*) -(L^R) 

Im(y4oo^i*) + (L R) 


( 6 ) 


as given in Ref. 14 , where an implicit integration over the mass of the K 7r~^ system is 
assnmed. 


The fraction of the S-wave contribntion is dehned as: 

l aL \2 , \ aR |2 

l^ool + l^ool 


Fsid) = 


dr/dg2 


(7) 


Where dT/dg^ is defined as the total differential rate of both S and P-wave contributions, 
given by 


dP _ dPs dPp 

dg2 

= ^ |^L|2 ^ |^L|2 ^ |^L|2 + ^ Jiy 


( 8 ) 


3 Fitting for the K*^ amplitudes 

3.1 Infinitesimal symmetries of the angular distribution 

Ignoring the S-wave terms, the angular distribution of —)■ decays can be 

described by eleven angular observables (Jj) for each flavour. These observables 
are made up of bilinear combinations of the K*^ spin amplitudes and represent the 
“experimental” degrees of freedom. If the J* terms are all independent, the experimental 
degrees of freedom should match the number of amplitude components which represent the 
“theoretical” degrees of freedom. However, there are continuous symmetry transformations 
of the amplitudes that leave the decay rate invariant [^. In order for the degrees of 
freedom to match it is required that 

rij - rid = 2na - ris, (9) 
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where Uj is the number of Jj terms, the number of the relations between the Ji, Ua is 
the number of complex amplitude components and Ug is a number of continuous symmetry 
transformations of the amplitudes that leave the decay rate invariant. In the massless 
limit ^ 4m^), and ignoring scalar contributions to the dimuon system, there are four 
continuous symmetry transformations of the amplitudes {ug = 4) that leave each of the 


Ji, and therefore the decay rate invariant, (see |11,34 for a detailed discussion). Given 
Uj = 11 and Ha = 12, as specihed by Eq. there are three relations between the various 
Ji, yielding eight independent angular observables. 

The following continuous transformations of the amplitudes, leave the angular distribu¬ 


tion unchanged 11,34 


0^4* L 

0 




cos 9 — sin 9 
sin 9 cos 9 


cosh iu — sinh iu 
- sinh ioj cosh ioj 


rii 


where the basis vectors n* are dehned as. 


ny 




n± = 


Ai 

_AR* 


no = 




( 10 ) 


( 11 ) 


The components (pi and (pn are phase-rotations of the left- and right-handed amplitudes 
separately. The second and third matrices act as a transformation between the left- and 
right-handed amplitudes. 

The angular distribution is degenerate under these transformations of the amplitudes. 
A likelihood function including all twelve real amplitude components, would therefore 
exhibit a 4D hypersurface of continuous maxima in amplitude space, rendering useless 
the minimisation techniques for the determination of the amplitude components. The 
symmetries of the angular distribution allow for the transformation of the amplitudes to 
a particular basis, where four of the amplitude components are hxed to some arbitrary 
value at every point in g^. The choice of the basis, referred to as “basis-hxing”, lifts 
the degeneracy. For the basis-hxing to allow for a subsequent ht of the remaining eight 
amplitudes as a function of g^, it is required that the values for 0^, (j)^, 9, and oj exist for 
every point in g^; and the amplitudes in this transformed basis are slowly varying in g^, 
such that they can be described by a simple functional form. This second requirement 
restricts the g^ range where the amplitudes can be extracted. The presence of potential 
light resonances below ~ 1 GeVyc"^ and of cc resonances above 8 GeVyc'^ motivates the 
use of the resonance-free and theoretically preferred region of 1 < g2 < 6 GeVVch 

A previous study described in Ref. used the following basis-hxing 

Re(Af) = Im(Af) = Im(Af) = Im(Af) = 0. (12) 

This basis suffers from a rapidly varying behaviour in Im(Ag) at g^ ~ 2 GeV^c^, as shown 
in Fig. In Ref. 0. the problems caused by this discontinuity were avoided by ignoring 
the g^ region below 2.5 GeV^c^. 

Ignoring the g^ region below 2.5GeVyc^ is clearly highly undesirable. For the method 
that we propose here, we instead use, 

Re(A^) = Im(A^) = Im(AQ) = Im(A;f) = 0. (13) 
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Figure 1: Distribution of the SM Im(Ag) amplitude in the untransformed (solid black) and in 
the fixed-basis proposed in Ref. (dashed blue). The untransformed amplitudes are given by 
the EDS program j^, up to an overall scaling factor. 


The amplitudes in the improved hxed-basis then exhibit a slow varying behaviour in 
both in the SM, shown in Fig. [^as well as in a range of new physics models. The EOS 
program is nsed to generate the amplitudes in the original basis. 


3.2 Exact discrete symmetries 

In addition to the continuous transformations of the amplitudes that leave the angular dis¬ 
tribution invariant, the angular distribution is also invariant under discrete transformations 
of the amplitudes. Even after the basis-hxing, which rednces the number of amplitudes 
to eight, there is a discrete symmetry comprising a simnltaneons shift Ai —>■ —Ai for all 
b that leaves the angnlar distribntion invariant. This symmetry can be seen simply by 
inspecting Eq. and noting that even after the conditions of Eq. IR have been applied, 
all angnlar observables are still constructed out of products of spin-amplitudes in the 
fixed-basis. 


3.3 Approximate discrete symmetries 


The limited amount of signal candidates available in the experimental data, can give rise 
to approximate symmetries nnder discrete transformations of the amplitudes. The exact 
form of these approximate symmetries can depend on the basis-hxing transformations 


discussed in Sec. 3.1 Given the basis-hxing condition of Eq. (13), a clear example occurs 


in the transformed basis with the SM amplitudes, where the lack of right-handed currents 
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Figure 2: Distribution of the transformed SM K*^(892) spin amplitudes in the fixed-basis 
(dashed black line) which exhibit a smooth behaviour in The solid red line denotes the result 
of the fit of the dependent ansatz discussed in the Sec. 3^ The untransformed amplitudes 
are given by the EOS program [^, up to an overall scaling factor. Only the non-zero amplitude 
components in the fixed-basis are shown. 


can give rise to an approximate symmetry under the transformation 

Au -A^ 




Al 


(14) 
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The effect of this accidental approximate symmetry can be demonstrated by generating 
samples based on the SM and on a model with large right-handed Wilson coefficients. 
Figure]^ shows the effect of the discrete transformation of Eq. 14 on cosOi, both in the 
SM and the model with large right-handed cnrrents. The region where —vr/d < 0 < 7r/2 
and 0 < cos 9 k < 1 is considered, in order to reduce possible cancellation of terms 
arising from the integration of the angular distribution over 0 and cos 6 k- Figure shows 
that in the SM the angular distribution in cos 6^ is essentially indistinguishable under the 
transformation; whereas, in the model with right-handed currents, the angular distributions 
can be distinguished. It is therefore clear that the above transformation is an approximate 
discrete symmetry of the angular distribution only in the SM and in other models with no 
right-handed currents. 




Figure 3: Demonstration of the effect of the transformation of Eq. (14) in the SM (left) and in a 
model with large right-handed Wilson Coefficients (right). The amount of data corresponds to 
the expected number of ^ candidates in LHCb’s Run-11 dataset. The blue line 

denotes the model that the data is generated from. The red-dashed line denotes the model with 
the approximate symmetry transformation, as mentioned in the main text, applied. The angles 
0 and Ok are required to satisfy —7r/4 < 0 < tt/2 and 0 < cos Ok < 1- 


An additional approximate discrete symmetry exists under the transformation of the 
right handed amplitudes in the transformed basis 


A? ^ -A 


R 


( 16 ) 


As for the left-handed amplitudes, the transformation of Eq. (15) is an approximate 


discrete symmetry of the angular distribution only in the SM and in other models with no 
right-handed currents. 


3.4 Parameterised amplitudes 

In order to determine the dependent K*^ spin amplitudes, a parametrisation of the 
amplitudes in the fixed-basis needs to be employed. A three-parameter ansatz of the form 

A = a + I3q^ +'y/q^ (16) 
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for both the real and imaginary components of each left- and right- handed amplitudes is 
chosen. No attempt is made to interpret these a, (3 and 7 coefficients in terms of short- or 
long-distance parameters. The choice of this ansatz is justihed by htting the transformed 
spin amplitudes, as provided by the EOS program [^, in the SM and numerous other 
physics models, using the parametrisation described above. Figure shows the result of 
the SM ht. Any bias coming from this choice of ansatz will be much smaller than the 
statistical uncertainty of current and any foreseeable-future experimental measurements. 

The basis-hxing reduces the number of amplitude components that need to be deter¬ 
mined to eight per flavour. Considering that each such component is described by three 
parameters to account for the dependence, in total there are twenty-four amplitude 
parameters per flavour that need to be determined. This parameter counting ignores 
any S-wave amplitudes. Such amplitudes are discussed further in Sec. |3.5[ Alternatively, 


the model dependent assumption can be made, that the only weak phases present in 
the amplitudes come from the CKM matrix elements. This assumption leads to the 
and amplitudes being identical, as the diagrams with non-zero weak phases are 
Cabibbo suppressed. Accounting for the experimental angular convention of the decay rate 
described in Sec. the decay distribution of both the B^ and 5° decays can be described 
using a single set of amplitude parameters. The approaches with separate and identical 
B^ and B^ amplitude parameters, are both discussed below. 

3.5 S-wave contribution 

Previous studies have discussed both the potential size as well as the impact of the S-wave 
contribution in the angular analysis of 5°—)- decays [^[^|^. In particular, it 

has been shown that the number of signal candidates expected in LHCb’s Run-I dataset, 
ignoring the S-wave contribution can have a signihcant effect on some angular observables. 
It is therefore critical that the dependent S-wave amplitude components are also 
accounted for in the ht to the angular distribution of B^ —)■ decays. 

In this study, the S-wave amplitudes are included in the angular distribution of the 
signal based on Eq. [^and are treated as nuisance parameters in the ht. The K^n~ mass 
range considered corresponds to 100 MeV/c^ around the A*‘’(892) pole mass. The rriKiT 
dependence is accounted for by modifying each term in Eq. by 


Jij — 


gi{mK-K)g*i{mKn)dmK-n 


(17) 


where gi{mKn) represents the rnKn line-shape of either a P-wave or an S-wave amplitude. 
It is thus assumed that the amplitudes do not explicitly depend on rriKn- 

The dependence of the S-wave amplitudes used in the generation of the simulated 
events are calculated following Ref. [^. For simplicity, only the k( 600) is considered to 
contribute to the S-wave in the K^7t~ mass range considered. This means that for this 
analysis only the line-shape of the k is considered to contribute to the S-wave giimKir)- 
The line-shape of both the «:(600) and the A**’(892) are taken as relativistic Breit-Wigner 


distributions with mass and width parameters as given in Ref. 40 . The form-factor of the 
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Table 1: Summary of the various integrals of the S- and P- wave line shapes that are used 
both in this study. The integral is performed in the rnKn range [796, 996] GeV/c^. 


Term 

Value 

/ \gK*o\‘^dmKn 

0.80 

f Ig^l^dniK^ 

0.18 

J" QkQ 

0.22 - 0.23f 


k( 600) is taken from Ref. J^. The values of the corresponding f gi{mKTT)gj{rnKw)dmKTT 
terms are shown in Tab. [Tjin the SM, the resulting value of Fs as a function of is 
shown in Fig. Using this simplistic approach, the predicted value of Fs is similar to the 
values obtained using more sophisticated treatments, such as those of Refs. 


37,38 


Given the size of current data samples, and the fact that the S-wave fraction is 
expected to be small (0(10%)), the dependence of the S-wave amplitudes Rqo^, can be 
approximated to be constant as a function of when performing hts to the data. This is a 
good approximation since the g^ shape of the S-wave amplitudes is expected to be similar 
, which is approximately constant in the region 1 < g2 < 6 GeV7cU 


to that of 


39 



Figure 4: Estimate of the S-wave fraction Fs in the SM as a function of g^ using the simplistic 
approach described in the text. 
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Table 2; Summary of the signal and background yields used in the generation of the 
simulated data samples. Events are assumed to be equally distributed between and 
candidates, both for signal and background. 


Sample 

Yields 

Run-I Run-II 

Signal 

Background 

600 2400 

500 2000 


3.6 Determining the amplitudes 

The stability and sensitivity of a £t for the dependent amplitudes is determined using 
simulated data with sample sizes equivalent to those expected at LHCb during Run-I 
and Run-II of the LHC. Estimates of the signal and background yields are taken from 
Refs. 15 33 and scaled linearly. Table [^summarises the signal and background yields used 
in the generation of the simulated data samples. The angular distribution of the signal is 
described using Eq. (j^, where the dependent amplitudes are again calculated using the 


EOS program for both the SM and new physics models. The angular distribution of the 
background is both generated and described in the ht as a product of four one-dimensional 
functions, each describing the dependence to the three helicity angles and g^, as shown in 


Eq. (18), 


d^r[Bkg] 


d cos 9^ d cos 9k d0 dg^ 


= /(cos6 *£) X g{cos9K) x h(0) x /(g^) , 


(18) 


where /, g, h, I are first order polynomials. 

Eor each dataset, the amplitude coefficients of Eq. IT are determined using an extended 
maximum likelihood ht. The probability distribution functions for the signal and back¬ 
ground, Psig(Bkg), are formed from the decay rate functions of Eq. [^and Eq. 18 respectively. 
The signal amplitudes, which make up the various Jj factors of Eq. [^ are written in terms 
of the three-parameter ansatz of Eq. The likelihood function is thus 


^Dat 

- logT = ^ -log[A^sig(aj,/5i,7i)^Sig(cos6'£,cos6'x,(/>,g^) + 

i 

Affikg.PBkg(cos6'£,cos6'x,0,g^)] + 

-iVcat log[Nsig{aj,Pj, 7j) A^Bkg] + [Nsig{aj,Pj, + iVekg], (19) 

where A^oat is the total number of events in the dataset, iVekg is a parameter in the ht 
that gives the number of background events, and iVsig(aj,y^) is the number of signal 
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events written in terms of the integrated signal decay rate: 


f^ji 1j) 


N, 


Dat 


»1 /•! /-TT /•6GeV2/c'‘ 


d^r[Sig] 


dcosOgdcosOK d(j)dq‘ 


d cos 9( d cos 9k dcj) dq^, 

( 20 ) 

where Nuat represents the average expected number of signal and background events for 
a given LHCb data sample, and Ag^ denotes the range in consideration, in this case 
5 GeV7ch 


4 Results 

An ensemble of 10^ simulated data sets is generated containing signal and background 
events as described in Sec. 3^ A maximum likelihood ht is performed to each of the data 
sets, to extract the g^ dependent P- and S-wave spin-amplitudes. Therefore at a given 
value of g^, 10^ determinations of each amplitude and thus of each angular observable are 
performed. The results of the hts for the g^ dependent amplitudes are presented in two 
scenarios. Firstly for a sample size equivalent to the full LHC Run-II data sample expected 
to be collected by LHCb, where the amplitudes of the and are both extracted 
without any model dependent assumptions (Scenario-!!); and secondly, for a sample size 
equivalent to the data sample collected by LHCb during Run-I of the LHC, where it is 
assumed that all weak phases of the amplitudes can be safely neglected (Scenario-I). In 
the latter case, the model dependent choice allows the sensitivity to the CP-averaged 
observables. Si, to be maximised, with the smaller sample that will be available from 
Run-I of the LHC. 

The resulting g^ dependent P-wave P° amplitudes, obtained from fits to an ensemble of 
simulated data under Scenario-!!, are shown in Fig. At a given point in g^, the 68% and 
95% conhdence intervals can be computed. Connecting these points at different g^ values 
gives the statistical uncertainty on the amplitudes as a function of g^. A clear degeneracy 
is observed under reflections about the x-axis. This effect is a consequence of the discrete 
symmetry Ai —)■ —Ai, as discussed in Sec. 3.2 Civen the observable quantities are bilinear 


combinations of the amplitudes, there is no corresponding degeneracy in any observable. 

The complete set of CP-symmetric and CP-asymmetric observables, dehned in Sec. 
can be constructed out the B^ and P° amplitudes. A subset of these observables are 
shown in Figs. |^and[^ for Scenario-!! and Scenario-I respectively. The 68% and 95% 
bands are given as in Fig. 

In a minor part of the g^ range, a bias at the level of 0.5cr is apparent in some 
observables, most notably 5*4. The 5*4 observable is sensitive to the interference between 
the Ay' and Ag amplitudes which in the fixed-basis is given by Re(A|'AQ). This bias arises 
from the approximate discrete symmetry discussed in Sec. |3.3[ This accidental symmetry 
is more prone to occur in models where the right-handed Wilson Coefficients are zero, such 
as the SM. As this is an effect of the angular distribution, this type of bias would have to 
be taken into account for any fitting method employed, whether binned or unbinned in g^. 
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q2 (GeV^/c^) (GeV^/c^) 


Figure 5; Distributions of the amplitudes as a function of resulting from 10“^ fits to 
generated signal and background simulated data under Scenario-II. The green and yellow bands 
correspond to the 68% and 95% confidence intervals respectively. The dotted line denotes the SM 
prediction as given by the EOS program [^. The magenta solid line denotes the most likely value 
resulting from the ensemble of fits. A discrete symmetry of reflections about zero is observed 
due to the fact that the Ji terms are bilinear coefficients of the amplitudes. 


4.1 Uncertainty estimation 

The probability distribution function of the signal decay remains positive dehnite for all 
values of the amplitude coefficients. This fact, conpled with the expected LHCb sample 
sizes for Scenarios-I and II, means that the likelihood surface, constructed out of the data 
and the probability distributions of the signal and the background, should be a good 
estimator of the statistical uncertainty of the amplitnde coefficients. 
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Figure 6: Subset of observables constructed out of the helicity-amplitudes, resulting from fits to 
ensembles of simulated data, under Scenario-II. The meaning of the various bands and curves is 
given in Figj^ 


Table summarises the pull mean and width of the amplitudes in the hxed-basis for 
Scenario-I at a point in q^. In order to discern any statistical bias from the bias originating 
from the approximate symmetries discussed in Sec. 3^, a point in is chosen such that 
minimises the bias from the approximate symmetries. 

The pull means and widths are largely consistent with zero and unity, respectively, 
indicating that the likelihood is a good estimator of the uncertainty of the amplitudes. 
The residual bias in arises from the additional approximate symmetry between 
and A^ discussed in Sec. 3.3[ Figure]^ shows that there is no point in q^ where the bias 
from the approximate symmetry can be removed for all the amplitudes simultaneously. 

A subset of two dimensional profile-likelihood distributions is shown in Fig. The profile 
likelihood is obtained by scanning over the two parameters in question, and minimising 
the likelihood over the rest of the parameters at each point. These prohle likelihoods can 
be used to determine the conhdence regions of the amplitude coefficients and therefore 
of the observables. The error matrix of the ht is a good approximation to the likelihood 
surface at the level of around 15%. 

For a given dataset, a prediction of the spin-amplitudes as a function of q^ can be 
obtained for a particular value of the Wilson Coefficients and B —)■ K* form factors. In 
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Figure 7: Subset of observables constructed out of the helicity-amplitudes, resulting from fits to 
ensembles of simulated data, under Scenario-I. The meaning of the various bands and curves is 
given in Figj^ 


turn, each amplitude can be expressed in terms of the three parameter ansatz of Eq. (16), 


Table 3: Means and widths of the pull distributions of the P-wave amplitudes at = 2.4 GeV^/c^, 
obtained from fits to ensembles of simulated data samples in Scenario-I. The pull is defined as 
(Fit-SM)/(TMeas. where UMeas. is the error on the measured quantity obtained using the error 
matrix of the fit. The deviation of the pull mean of Re(M|j^) from zero arises due to the residual 
bias from the approximate symmetry of the angular distribution. 


Parameter Pull mean Pull width 


Re{A^)B^ 
lm{A^)B^ 
Re{Af)B^ 
lm{Af)B^ 
Re{A^)B^ 
lm{A^)B° 
Re (71^)50 


-0.03 

± 

0.02 

0.00 

± 

0.02 

0.01 

± 

0.02 

0.22 

± 

0.02 

-0.02 

± 

0.02 

0.02 

± 

0.02 

-0.04 

± 

0.02 

-0.05 

± 

0.02 


0.97 ±0.03 
1.01 ±0.03 
1.02 ±0.03 
0.90 ±0.03 
0.94 ±0.03 
0.97 ±0.03 
0.95 ±0.03 
0.94 ±0.03 
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Figure 8: Subset of two dimensional profile likelihood distributions from a fit to a single simulated 
dataset in Scenario-I. The red ellipses denote the error and the correlation obtained from the 
error matrix of the ht. 


by transforming to the fixed-basis, nsing the procedure detailed in Appendix]^ By fitting 
each transformed amplitude with the ansatz of Eq. (16), the coefficients a, (3 and 7 can 
be determined. The best fit point and the error matrix of the fit to the data, can then be 
used to place constraints on the Wilson Coefficients using a procedure like that employed 


in Refs. 
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5 Sensitivity to new physics 


The expected sensitivity to the effects of new physics (NP), neglecting theory uncertainties, 
is estimated by generating a large number of simulated data samples according to the 
SM (null dataset) and a NP model where 6Cq = —1.5 (test dataset). The choice of 
the NP model is motivated by recent results from global fits to 6 —)■ si£ and b ^ sj 
measurements [6 17-28 which are dominated by LHCb’s anomalous results in the angular 
distribution of 5°—decays 15,16. The EOS program is used to generate 
simulated data from these two models using their central value predictions. Two fits are 
then performed to an ensemble of simulated datasets generated with the SM and the 
NP model. In the first fit, the amplitude parameters are fixed to their SM values (null 
hypothesis). In the second fit, the amplitude parameters are fixed to the values given by 
the model with 6Cg = —1.5 (test hypothesis). The background components and yields are 
treated as nuisance parameters and are left floating in each fit. The S-wave contribution 
in the Kn system of —)■ decays is less understood theoretically than the 

dominant P-wave part. In order to correctly model the S-wave component of the Kn 
system, experimental input is required. Therefore, for these sensitivity studies, the S-wave 
amplitudes are treated as nuisance parameters in the fit. 

The test statistics are defined as 


Q 

Q 


SM 

NP 


2{NLLif, - NLLf^lf,) 
2(NLL,^i^-NLL(l^,), 


( 21 ) 


where ^^^nuiue^t corresponds to the negative log likelihood value of the null or test 
hypothesis on a SM or NP simulated dataset. The expected sensitivity to a model with 
SCg = —1.5 is then estimated by counting the fraction of the toy simulations with a value 
of where is the median of the distribution. Figure shows the 

distribution of the test statistic for both the SM and NP simulated data samples in fits to 

and B^ candidates separately, using a sample size equivalent to that expected at LHCb 
during Run-I of the LHC. The probability for the SM sample to fluctuate such that it 
gives a test statistic as low or lower than the median of the NP sample (i.e 
corresponds to a significance of 6.5a. This significance is obtained for an idealised model 
of the experimental data and does not account for the theoretical uncertainties associated 
with the translation from Wilson coefficients to amplitudes. 

For comparison with methods used previously, the same procedure for estimating the 
expected significance can be performed for fits directly to the CP-averaged observables. 
Si in bins of For these fits, three bins between 1 < g^ < 6 GeV^C are chosen as 
(1, 2.7), (2.7,4.3), (4.3, 6) GeV^C. In this binned approach a combined significance of 5.0(j 
is obtained. Therefore, fitting for the g^ dependent K*^ spin-amplitudes, separately for 
the B^ and the P°, results in a 30% improvement in the expected sensitivity. Equivalently, 
to get the same sensitivity as the amplitude method, a binned fit to the CP-averaged 
observables would require an additional 70% of integrated luminosity. 

The inclusion of the Ktt system in an S-wave configuration introduces six additional 
observables, as shown in Eq. ([^, in contrast to four additional S-wave amplitude compo- 
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nents. Approximately half of the increase in sensitivity detailed above, can be attributed 
to the reduced number of S-wave-related nuisance parameters present in the amplitude 
fits, and the other half owes to the intrinsic sensitivity of the method. 



2(NLL^^^^-NLL„^„) 


Figure 9: Distribution of the test statistic ]\gp points) and SM 

(blue open points) simulated data samples. The S-wave components are treated as nuisance 
parameters and are therefore floating in the fit. The dashed red line denotes the median of the 
test statistic from fits to the NP sample. The solid blue curve is a fit of a Gaussian distribution 
to the distribution of the test statistic from fits to the SM sample. 


6 Conclusions 

In summary, a method of analysing —)■ decays is presented that allows the 

determination of all of the K*^ amplitudes as a parametric function of q^. The method is 
applicable with the data sample that is already available at the LHCb experiment and 
works in the region where the amplitudes can be described by a simple functional form, 
f < < 6 GeVyc"^. The bias coming from the choice of parameterisation is much 

smaller than the statistical uncertainty of current and any forseeable-future experimental 
measurements. 

The method overcomes several shortcomings of previous methods: As all the K*^ 
helicity amplitudes are determined from a single ht to data, the full correlations between 
experimentally determined quantities can be obtained, improving the sensitivity to new 
physics; fitting for the amplitudes enables all of the symmetries of the angular distribution 
for a Ktt system in a P- and S-wave state to be accounted for, giving increased experimental 
precision compared to previous approaches which retain redundant degrees of freedom; 
the method avoids the integration of theoretical predictions over experimental q^ bins. 
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enabling the fnll exploitation of the cancellation of the form factors at leading order and 
therefore fnrther improving the sensitivity to new physics. 

The continuous symmetry transformations that can be applied to the amplitudes enable 
a modihed basis to be specihed, in which a number of amplitude components are set 
to zero. The basis considered by previous studies 11 gave a discontinuous shape in 


rendering the region < 2.5 GeV^c^ unusable. Here, a new choice of basis is presented, 
which leaves the amplitudes smoothly varying in the entire range 1 < g^ < 6 GeV^c"^ for 
the SM and for a wide range of new physics models. The present analysis also highlights 
new approximate discrete symmetries that are manifest for data samples on the order of 
that collected at LHGb during Run-I. These approximate symmetries are more prone to 
occur in models like the SM, where there are no right-handed currents. 

In order to illustrate the options that are available to ht a given dataset, the sensitivity 
of the method is presented in two different scenarios. For a sample equivalent to the LHGb 
Run-I dataset, a model-dependent assumption is made that the only weak phases present 
in the amplitudes come from the GKM matrix elements. Given diagrams with non-zero 
weak phases are Cabibbo suppressed, this assumption results in and amplitudes 
which are identical. The CP-averaged observables can then be determined with greater 
precision than would be possible if instead both the CP-averaged and CP-asymmetric 
observables were determined. For a sample equivalent to the LHGb Run-H dataset, the 
results are determined without this assumption and the sensitivity to both CP-averaged 
and CP-asymmetric observables is presented. 

The recent anomalous result in the angular distribution of P° motivates 

the consideration of a new physics scenario with 6Cq = —1.5 15 . In such a scenario, the 


new method presented is 30% more sensitive than a three g^-bin fit to the CP-averaged, Si, 
in the region 1 < g^ < 6 GeV^C, improving the discrimination between this scenario and 
the SM from 5.0 to 6.5 a. An additional ~70% of integrated luminosity would therefore be 
required in order for a £t to the CP-averaged observables to achieve the same sensitivity 
as the fit to the K*^ spin-amplitudes. This improvement arises from the treatment of the 
g^ dependence of the amplitudes as a continuous function, and the use of all symmetry 
relations of the angular distribution to reduce the number of independent parameters that 
describe the P° —)■ decay. 
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Appendix 


A Amplitude transformation and parametrisation 


The spin-amplitudes in the fixed-basis can be determined by applying Eq. ( |10[ ) to a set of 
amplitudes in the original basis. The transformation angles, 0/,, 0/^, 0, a; can be determined 
in terms of the original amplitudes, by requiring that in the fixed-basis Eq. (13) holds, 
and solving a system of four non-linear equations. The analytical expressions of the 
transformation angles are given in Ref. and shown below 


tan 2u 
tan 6 ' 

tan 
tan (pfi 


Im(/l«)Re(/l^) + (L ^ fi) 

\A§\^ - 

Re(R^) -I- Im(74o) tana; 

—Re(74Q) -|- Im(74^) tana; 

^(Aq) -|- Im74^tan6' — [Re(74^) — Re{AQ) tan^] tana; 

—Re(74Q) -|- Re(74^) tan^ -|- [Im(74^) -|- Im(74Q) tan^] tana; 
Im(R;5^) -|- Im(74^) tan 6 — [Re(74;5l) ~ Re(74;5^) tan 6] tan u 
—Re(74;5^) -I- Re(74^) tan6' -|- [Im(R^) -|- Im(74;5^) tan^] tana;’ 


where are the K*^ amplitudes in the original basis. 

A C++ library will be provided shortly, which will transform a given set of amplitudes 
to the fixed-basis, and perform a fit in to determine the coefficients a, f3, 7 of the 
ansatz of Eq. (16). In addition, provided there is a consensus in the experimental 
community, the likelihood function of Eq. (19) could be made publicly available, along 
with tools that minimise the likelihood over the nuisance parameters, for a given set of 
amplitudes. 
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